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A HYBRIDIZED FORMULATION FOR THE WEAK GALERKIN 
MIXED FINITE ELEMENT METHOD 

LIN MU*, JUMPING WANGt, AND XIU YE* 


Abstract. This paper presents a hybridized formulation for the weak Galerkin mixed finite 
element method (WG-MFEM) which was introduced and analyzed in m for second order elliptic 
equations. The WG-MFEM method was designed by using discontinuous piecewise polynomials on 
finite element partitions consisting of polygonal or polyhedral elements of arbitrary shape. The key 
to WG-MFEM is the use of a discrete weak divergence operator which is defined and computed by 
solving inexpensive problems locally on each element. The hybridized formulation of this paper leads 
to a significantly reduced system of linear equations involving only the unknowns arising from the 
Lagrange multiplier in hybridization. Optimal-order error estimates are derived for the hybridized 
WG-MFEM approximations. Some numerical results are reported to confirm the theory and a 
superconvergence for the Lagrange multiplier. 

Key words, weak Galerkin, finite element methods, discrete weak divergence, second-order 
elliptic problems, hybridized mixed finite element methods 

AMS subject classifications. Primary, 65N30, 65N15; Secondary, 35J20, 76S05, 35J46. 


1. Introduction. In this paper, we are concerned with new developments of 
weak Galerkin finite element methods for partial differential equations. Weak Galerkin 
(WG) [211 UHl EQI HH] is a generic finite element method for partial differential equa¬ 
tions where the differential operators (e.g., gradient, divergence, curl, Laplacian etc) 
in the variational form are approximated by weak forms as generalized distributions. 
This process often involves the solution of inexpensive problems defined locally on 
each element. The solution from the local problems can be regarded as a reconstruc¬ 
tion of the corresponding differential operators. The fundamental difference between 
the weak Galerkin finite element method and other existing methods is the use of 
weak functions and weak derivatives (i.e., locally reconstructed differential opera¬ 
tors) in the design of numerical schemes based on existing variational forms for the 
underlying PDE problems. Weak Galerkin is, therefore, a natural extension of the 
conforming Galerkin finite element method. Due to its great structural flexibility, the 
weak Galerkin finite element method fits well to most partial differential equations by 
providing the needed stability and accuracy in approximation. 

The goal of this paper is to develop a new computational method which reduces 
the computational complexity for the weak Galerkin mixed finite element methods 
[19| by using the well-known hybridization technique [141 [1]. Our model problem 
seeks a vector-valued function q = q(x), also known as flux function, and a scalar 
function u = u(x) defined in an open bounded polygonal or polyhedral domain 12 C 
(d = 2, 3) satisfying 

(1.1) aq-|-Vu = 0, V • q =/, in D 
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and the following Dirichlet boundary condition 


( 1 . 2 ) 


u = g, on dfl, 


where a = {aij{x))dxd € is a symmetric, uniformly positive definite 

matrix on the domain ft. A weak formulation for (1.1)-(1.2) seeks q G H(dw,fl) and 
u e L^(fi) such that 


(1.3) (aq,v) - (V • V, u) = -(g,v • n)an, Vv G H(div,n) 

(1.4) (V • q,u;) = (/,«;), Vw G L^(n). 

Here L^(fl) is the standard space of square integrable functions on fl, V • v is the 
divergence of vector-valued functions v on Q, H{div, fl) is the Sobolev space consisting 
of vector-valued functions v such that v G and V • v G i^(0), (•, •) stands 

for the L^-inner product in L^(n), and (•,-)an is the inner product in 

Conforming Galerkin finite element methods based on the weak formulation 
) and finite dimensional subspaces of H{div, 11) x (12) with piecewise polynomials 

known as mixed finite element methods (MFEM) [151 114] . MFEMs for (dH)- 
) treat q and u as independent unknown functions and are capable of providing 
accurate approximations for both unknowns [maiHllTllgllinillllllSllle]. ah the 
existing MFEMs in literature possess local mass conservation that makes MEEM a 
competitive numerical technique in many applications such as flow of fluid in porous 
media including oil reservoir and groundwater contamination simulation. 

Based on the weak formulation (|1.4[ ), a weak Galerkin mixed finite element 
method (WG-MFEM) was developed in [19] which provides accurate numerical ap¬ 
proximations for both the scalar and the flux variables. Like the existing MEEMs, the 
WG-MEEM scheme conserves mass locally on each element. But unlike the existing 
MFEMs, the WG-MFEM allows the use of finite element partitions consisting of poly¬ 
gons {d = 2) or polyhedra {d = 3) of arbitrary shape that satisfy the shape-regularity 
assumption specified in |19| . It should be pointed out that some similar features are 
shared by a number of recently developed numerical methods such as Virtual Element 
Methods mm, hybridizable discontinuous Galerkin methods m, and mimetic finite 
differences 011]. 

While weak functions and weak derivatives provide a great deal of flexibility for 
WG methods, they also introduce more degrees of freedom than the standard finite 
element method. The purpose of this paper is to develop a hybridized formulation 
for the weak Galerkin mixed finite element method m that shall reduce the compu¬ 
tational complexity significantly by solving a linear system involving a small number 
of unknowns arising from an auxiliary function called Lagrange multiplier. The La¬ 
grange multiplier is defined only on the element boundaries (also called wired-basket 
of the finite element partition). As a result, the linear system that costs the majority 
of the computing time depends on the dimension of the finite element space defined 
on the wired-basket. Optimal-order error estimates for the hybridized WG-MFEM 
(HWG-MFEM) approximations are established in several discrete norms. Some nu¬ 
merical results are presented to demonstrate the efficiency and power of the hybridized 
WG-EEM. 

Throughout the paper, we will follow the usual notation for Sobolev spaces and 
norms [niaiTj. For any open bounded domain D C with Lipschitz continuous 
boundary, we use || • ||s,d and | • |s,d to denote the norm and seminorms in the Sobolev 
space for any s > 0, respectively. The inner product in H^{D) is denoted by 


are 
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(•, ■)s,D- The space H^{D) coincides with L'^{D), for which the norm and the inner 
product are denoted by || • ||d and respectively. When D = Q, we shall drop 

the subscript D in the norm and inner product notation. 

The paper is organized as follows. In Section 2, we review the discrete weak 
divergence operator. In Section 3, we present a hybridized weak Galerkin mixed finite 
element method. Section 4 is devoted to a discussion of the relation between the 
WG-MFEM and its hybridized version. In Section 5, we derive an error estimate for 
the hybridized WG-MFEM. Finally, in Section 6, we report some numerical results 
that demonstrate the efficiency and accuracy of the hybridized WG-MFEM, including 
a superconvergence for the Lagrange multiplier. 

2. Weak Divergence. Let iF be a polygonal or polyhedral domain. A weak 
vector-valued function on K refers to a vector field v = {vo,vt,} with Vq S [L‘^{K)Y 
and Vb G \L?'{dK)Y . The first component Vq carries the information of v in K, and 
Vf, represents partial or full information of v on dK. The choice of the boundary in¬ 
formation that Vfc represents is problem-dependent. Note that Vf, may not necessarily 
be related to the trace of Vq on dK should a trace be well-defined. 

Denote by V{K) the space of all weak vector-valued functions on AT; i.e., 

(2.1) V{K) = {v = {vo, Vb} : Vq G [L\K)]\ v, G [L\K)Y}. 

A weak divergence can be taken for any vector field in V(K) by following the definition 
introduced in |19| . which we summarize as follows. 

Definition 2.1. For any v G V{K), the weak divergence of v is defined as 
a linear functional, denoted by • v, on F[^{K) whose action on each ip G H^{K) 
is given by 

(2.2) (Vu, • V, lp)k ■= -(vo, Vif>)K + (vb • n, tp)dK, 

where (•, ■)k and (•, ■)qk stands for the Lfi-inner product in L‘^{K) and L‘^{dK), re¬ 
spectively. 

The Sobolev space [H^ {K)Y can be embedded into the space V(K) by an inclusion 
map iy : [F[^{K)Y -G V{K) defined as follows 

w(q) = 

With the help of the inclusion map iy, the Sobolev space [H^{K)Y can be viewed 
as a subspace of V{K) by identifying each q G [H^{K)Y with *y(q)- Analogously, a 
weak function v = {vo,Vb} G V{K) is said to be in [H^{K)Y if it can be identified 
with a function q G \H^{K)Y through the above inclusion map. It is not hard to see 
that Vuj • V = V • V if V is a smooth function in [F[^{K)Y. 

Next, we introduce a discrete weak divergence operator by approximating (V^, • ) 
in a polynomial subspace of the dual of F[^{K). To this end, for any non-negative 
integer r > 0, denote by Pr{K) the set of polynomials on K with degree r or less. 

Definition 2.2. A discrete weak divergence (Vw.r') is defined as the unique 
polynomial (S/w,r ’ v) G Pr{K) satisfying the following equation 


(2.3) 


(V^,,r ■ V, (j))K = -(vo, \7(j))K -f (Vb • n, (tfidK, V(/) G Pr{K). 
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3. Hybridized WG-MFEM. Let 7^ be a finite element partition of the domain 
n consisting of polygons in two dimensions or polyhedra in three dimensions satisfying 
the shape-regularity condition specified in m- For T G Th, denote by hx its diameter 
and h = maxTeTh the meshsize of 77- The set of all edges or flat faces in Th is 
denoted as £h, with the subset = £h\d£l consisting of all the interior edges or flat 
faces. 

Let fc > 0 be any integer. On each element T gTh, we denote by n the outward 
normal direction on the boundary dT, and define two local finite element spaces 
Wk+iiT) = Pk+iiT) and 

(3.1) Vfc(T) = {v = {vo,Vb} : Vq S [Pk(T)f, Vb\e = Vbn, Vb G Pk{e), e C dT}. 

On the wired-basket £h^ we introduce a finite element space using piecewise polyno¬ 
mials of degree k: 

(3.2) A;, = {A : Ale e Pfe(e), Ve e £h}. 

Let A° C Afi be the subspace consisting of functions with zero boundary value 

(3.3) A° = {A e A;,; A|e = 0, Ve C 50}. 

We further introduce an element-wise stabilizer as follows 


(3.4) 


ST(r, v) = hrUro - n) ■ n, (vq - v^) • n)aT, 


r, V 


e VkiT). 


Denote by the discrete weak divergence operator in the space 14 (T) 

computed by using (2.3) on each element T. For simplicity of notation, from now 
on we shall drop the subscript fc -I- 1 from the notation for the discrete weak 

divergence. 


Hybridized WG-MFEM Scheme 3.1. For an approximate solution of (1.1)- 
find q/i = {qo,q6} G Vk{T), Uh G Wk+i{T), Xh G Ah such that Xh = QbP on 
d£l and 


(3.5) ST(q/i, v) + (aqo, vo)t - (v^, • v,u/,)t =-(A/i, Vb • n)aT, Vv e 14(T), 

(3.6) (Vu, • q/i,-u;)T = (/, w)t, VwGWk+iiT), 

(3.7) (qb • n, <l))aT = 0, G A^. 

TdTh 


Here {Qbg)\e S i4(e) is the Lf projection of the boundary value u = g on the edge 
e C d£l. 


The approximation for the primal variable u is given by Xh on £h and Uh on 
each element T. The flux q is approximated by qo on each element, and qb is an 
approximation of the normal component of q on the element boundary. 


The hybridized WG mixed finite element scheme (3.51-(3.7) is an analogy of the 
hybridized mixed finite element method [14] (see mtiMi for more details). Note 
that the approximation q^ and Uh are defined locally on each element T G Th so 
that qb assumes multi-values on each interior edge. More specifically, unlike the WG 
mixed finite element method introduced in [T5], the hybridized WG-MFEM scheme 
(3.5)-(3.7l does not assume the “continuity” (or single-value) of qb on each interior 
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edge e € The new variable Xh, known as the Lagrange multiplier, wa s ad ded here 
to provide the necessary “continuity” of qb on through the equation (3.7). 


In the rest of this section, we shall reformulate the problem (3.5)-(3.7) by elimi¬ 


nating the unknowns q/j and u/,. The resulting linear system is symmetric and positive 
dehnite in terms of the variable Xh- Readers are referred to [T] for a similar result for 
the standard mixed finite element method. 

Denote by X4h the collection of all the local hnite element spaces: 


Mh= (g) Vk(T) X Wk+i(T). 

Ten 

For any given 9 S Ah, let {q®,M®} S Mh be the unique solution of the following 
problem 

(3.8) ST(q®,v)(aq^,vo)T - (V^ • v,m^)t =-(6»,Vb • n)aT, Vv S I4(T), 

(3.9) (V^-q^rc)T = 0, yw€Wk+i{T). 

By setting 

H(0) :={q^,q^,<} 

we see that El defines a linear operator from Ahto Aih- For convenience, we decompose 
the operator M into three components El = {Hq,E tf,, Hu} so that 

q^=Ho(0), q^=Hb(0), 4 =Hu(0). 


The first two components of H can be grouped together to give the following operator 

Hq = {Ho,Hb}. 

Note that the third component Hu(6*) is a piecewise polynomial of degree fc -I- 1 which 
is indeed a discrete harmonic extension of 9 in the domain D from the wired-basket 
Sh- 

With the help of the linear operator H, we can reformulate the hybridized WG- 
MFEM scheme (3.5)-(3.7) into one involving only the unknown variable on the wired- 
basket £h- 

Theorem 3.2. The solution Xh € Ah arising from the hybridized WG-MFEM 
scheme (3. 5)- (3.1) is a solution of the following problem: Find 9 G Ah satisfying 
Xh = Qbg on dXl and the following eguation 

(3.10) ^ {sT(Hq(A,,),Hq((^)) + (aHo(A„),Ho(<(.))T} = (/,Hu(0)), V<^ e A^. 

Ten 


Proof. Note that the equation (3.7) has test function (f G A°. Owing to the 
operator H, by letting 9 = cf in (|3.8|), we arrive at 


-(</), vb • n)aT = ST(Hq(()i), v) -h (aHo((/i), vo)t - (V^, • v,H„((/)))t 


for all V G I4(r). In particular, by setting v = q/, (i.e., the solution of (3.5)-(3.6)), 
we obtain 


-{(j), qb • n)oT = ST(Hq(0),q/,) -h (aHo((()), qo)r - (V^, • qb,H„((^))r 
= ST(q/i,Elq((^)) -b (aqo,Ho((()))T - (/,H„(0 ))t, 
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where we have used (3.6) in the second line. By summing over all the element T we 
have from (13.7D that 


(3.11) 


^ {sTicih,n^m + (aqo,Ho(.^))T} = V<^ e A^. 

TGTh 


Next, by using (3.9) and then (3.5) with v = we obtain 

ST(q/i,IHIq((^)) + (aqo,IHIo((/)))T 
= ST(qh,IHIq(^)) + (aqo,IHIo((/)))T - (V^, • B.q{(l)),Uh)T 
= — (A/i, • n)gr 

= ST(IHIq(A/i),IHIq((/))) + (aIHIo(A/i),IHIo(()>))r — (V^, • EIq((()),EI„(A;i))T 

= ST(Hq(A;,),Hq((/.)) + (aHo(A„),Ho(</.))T. 

Substituting the above equation into ( |3.11 ) yields 

^ {sT(Hq(A,),Hq(,^)) + (aHo(A^),Ho(<).))T} = ^4> G K. 

TGTh 


which is precisely the equation (3.10). □ 


Concerning the solution existence and uniqueness for the reduced system (3.10), 
we have the following result. 

Theorem 3.3. There exists one and only one Xh G A;, satisfying the equation 
( 3l^ and the boundary condition Xh = Qtg on dfl. 

Proof. Since the number of unknowns equals the number of equations, it is suf¬ 
ficient to verify the solution uniqueness. To this end, let A^*^ G Ah be two solutions 

I I 

of (3.10) satisfying the boundary condition Xf^ = Qbg, * = 1,2. It is clear that their 
difference, denoted by cr = A^^^^ — Ajf\ is a function in A° and satisfies 

^ {sT(Hq(a), Hq(</>)) + (aHo(a), Ho(<)>))t} = 0, V</) G A°. 

TdTh 

By setting (f = a we arrive at 

ST(IHIq(CT),IHIq(cr)) -b (aHo( ct), Hq( cr ))t = 0, V T G Tft. 

It follows that = 0 on each element T and • n = Ho (cr) • n on dT, which 

implies H{,(cr) = 0 on dT. Now using (3.8) we have 


(3.12) 


(Vuj • v,H„(cr))T = (cr, Vb • n)aT, Vv G Vfc(r), T G %. 


The definition of the discrete weak divergence can be applied to the left-hand side to 
yield 

(3.13) - (vo, VH„(cr))T-b (vb • n,H„(cr))aT = (cr,V{, • n)OT, Vv G Vfc(r). 

Thus, by setting Vb = 0 and varying Vq we obtain 

VHi,(cr) = 0 ==b Hu(cr) = const, in T G 7b, 


which, together with varying Vb in ( |3.13[ ), leads to 
that both 


j(cr) = cr on dT. It follows 


(cr) and cr are constants in the domain and vanishes on the boundary. 
Hence, a = Hii(cr) = 0. This completes the proof of the theorem. □ 
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4. Equivalence between HWG-MFEM and WG-MFEM. The goal of this 
section is to show that the hybridized WG-MFEM approximate solution arising from 
) coincides with the solution from the corresponding WG-MFEM scheme in¬ 
troduced as in m- To this end, consider the finite element space 14 = ^reTh 
For any interior edge/face e € denote by Ti and T 2 the two elements that share e 
as a common side. Recall that functions v = {vo,V{,} € 14 have two-sided values on 
e for the component v;,: one comes from the element Ti and other from T 2 . Define 
the jump of v along e G as follows 

(4-1) [v]e = ■ ni + ^b\aT^ ■ n2, 

where is the outward normal direction on e as seen from the element Ti, i = 1,2. 

Let fc > 0 be any integer. Following m, we introduce two weak finite element 
spaces 

Vh = {v = {vo,Vb} : V G 14, |v]e = 0, Ve G S^}, 

Wh= ^ Wk+i{T). 

TGTh 

On each element T £ Th, introduce three bilinear forms 


(4.2) 

arir, v) = (aro, vo)t, 

(4.3) 

&t(v, w) = (Vu, • V, w)t, 

(4.4) 

ct(v, a) = (vf, • n, a)Qj,, 


for V = {vo,Vf,},r = {ro,rf,} G Vk{T), w G Wk+i{T) and a G A/j. In addition, set 

as^rir, v) = arir, v) -b srir, v). 


(3.5)-(3.7 


WG-MFEM Scheme 4.1. For an approximation of the solution of (1.1)- 
(fl.ll), find = {qo,qfc} € Vh and Uh G 114 such that 


(4.5) 

(4.6) 


X! “s.T(q/i,v) - X! bT{'v,Uh) = -{g, Vh • n)9s.j, Vv G Vh 

TGTfc TGTh 

^ briqh^w) = {f, w), Vic G IF/,. 
Tan 


The following Theorem shows that the WG-MFEM is equivalent to the hybridized 
WG-MFEM. 

Theorem 4.2. Let {qh,Uh, G A 4 h x A/, be the hybridized WG-MFEM ap¬ 
proximation of (1.1)-(1.2) obtained from ^3.1 


___ ,__ __ and (qhjUh) G Vh X Wh be the 

WG-MFEM approximation of O - (7.^ arising from Then we have 

qh = qh o,Td Uh = Uh- Moreover, the hybridized WG-MFEM scheme (3.5)-(3.1) has 
one and only one solution. 

Proof As sume that (q/,, it/,. A/,) G Mh x A/, is a solution of (3.5)-(3.7). From the 
equation (3.7) we see that [q/,| = 0 on each interior edge e G S'f. Hence, q/, G F/,. 
Moreover, for any v G F/, so that |v/,[|e = 0 with e G 5°, we have 


(A/„v 6 • n)aT = {g,^b ■ n)afj. 

Ten 


(4.7) 


















Therefore, for v G Vh and w G Wh, the equations (|3.5|)-(3.6| leads to 


as.T(q/i,v)- Y bT{^,Uh) = -{g,Vb ■ n)g^, 

TGTh TeTh 

Y bTi(lh,w) = (/, w), 

Ten 

which shows that {qfi,Uh) is a solution for the WG-MFEM. The solution existence 
and uniqueness for the WG-MFEM [19] then implies q;i = q;i and Uh = Uh- 


To show that the hybridized WG-MFEM scheme (3.5)-(3.7) has a unique solution, 
assume that ^ ■^h x * = 1)2, are two solutions. Their difference, 

[cih,Uh,Xh) := - o^h\u^h - 


must satisfy ( |3.5| )-(3.7 1 with / = 0 and g = 0. From the first part of this theorem, we 
see that 


q/i = 0, Uh = 0. 


Thus, from (3.5) we have 


0 = (A/jjVb • n)QT, Vv G Vk(T), T gTh, 


which implies A/i = 0 on dT for any T G Th- This completes the proof of the theorem. 

□ 


The hybridized WG-MFEM scheme (3.5)-(3.7l consists of two parts: a local sys¬ 
tem (|3.5|)-(|3.6|) defined on each element T G Th and a global system (13.7|). This 


scheme has a reformulation of (3.10) which involves only the Lagrange multiplier Xh- 


It should be pointed out that the reduced system (3.10) is symmetric and positive def¬ 


inite. Since the Lagrange multiplier Xt is defined only on the element boundary as its 
unknowns, the size of the linear system for the hybridized WG-MFEM is significantly 
smaller than the one for the WG-MFEM scheme. Thus, the hybridized WG-MFEM 
is an efficient implementation of the WG-MFEM scheme proposed and analyzed in 

m- 


5. Error Estimates. 

MFEM scheme 
developed in [T 


Owing to the equivalence between the hybridized WG- 


(3.5)-(3.7|. 


3.5| )-(3.7) and the WG-MFEM scheme of [12], all the error estimates 
I can be applied to the approximate solution (q^,u/i) obtained from 
What remains to study is the convergence and error estimate for the 
Lagrange multiplier Xh- 

On each element T G Th, denote by Qq the projection from [L'^{T)Y‘ to 

[Pk{T)Y, by Qb the projection from L^(e) to Ffe(e), and by the projection 
from L'^(T) onto Pk+i{T). With the help of these operators, we define a projection 
operator Qh ■ V{T) -G VuiT) so that for any v = {vo,?;;,!!} G V{T) 


(5.1) 


Qh'v = {QoVq, {QbVb)n}. 


In the finite element space Wh, we introduce the following norm: 


\l,h 


Y 

T^Th 


+ h 


-^Y\ 

eeSh 


(5.2) 
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In the space 14, we use 


a*.T(v,v). 

TGTh 


For simplicity of notation, we shall use ^ to denote “less than or equal to” up to a 
constant independent of the mesh size, variables, or other parameters appearing in 
the inequalities of this section. 


Theorem 5.1. m Let { qh , Uh , Xh) G A4h x -Xh be the approximate solution of 
{ arising from the hybridized WG-MFEM scheme (3.5)-(3.1) of order k > 0. 

Assume that the exact solution (q,u) of (l.S)-(1.4) is regular such that u G 
and q S Then, we have 


(5.3) |||q/i - Qhq\\\ + \\uh - Qhu\\i,h < (||u|U +2 + |lq|U+i) ■ 


If, in addition, the problem 0 with homogeneous Dirichlet boundary condition u = 
0 has the usual -regularity, then the following optimal order error estimate in 
holds true: 


(5.4) 


||R/.-Q^R||</i'=+"(||R|U+2 + |lq|U+l). 


For any v = {vq, Vf,} S Vk{T) and w G H^{T), using the definition of Qh and the 
integration by parts we obtain 

(Vu, • V, <Qhw)T = -(vo, V(Q;iUi))t + (vf, • n, <Qhw)oT 

= (V • Vo, Qh'w)T + (Vb • n - Vo • n, Qhw)dT 
= (V • Vo, w)t + (vfo • n - Vo • n, <Qhw)dT 
= -(vo, \/w)t + (vo • n, w)oT + (vfc • n - Vo • n, Qhw)dT 
(5.5) = -(vo, \/w)t + (vo • n - Vf, • n, w- Qhw)dT + {^b ■ n, w)ot. 


By introducing the following notation 

iriw, v) = (vo • n - Vb • n, QhW - w)ot, 


we have from (5.5) that 
(5.6) 


(Vu, • V, Qhw)T = -(vo, Vw)t + (vfe • n, u;)aT + ^t(w, v). 


Lemma 5.2. Let {qh,Uh, Xh) G Aih x A/, be the approximate solution of {1. li¬ 
ft. 2) arising from the hybridized WG-MFEM scheme {3.5)-{3.7) of order k > 0. On 
each element T G Th, the following identity holds true 

ST(eh, v) + (aeo, Vo)r - (V^ • v, cuIt + {bh, Vf, • 

= sriQhW v) + ^t(m, v), Vv e Vfc(r), 

where 


— Qh^ q/i 5 bfi — Qb'^ Xf,, 

stand for the error between the hybridized WG-MFEM approximate solution and the 
L^ projection of the exact solution. 
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Proof. For any v = {vo,Vf,} G Vk{T), we test the first equation of (O) against 
Vq to obtain 

(aq,Vo)T + (Vw,Vo)r = 0. 

By applying the identity ( |5.6[ ) to the term (Vm,Vo)t we arrive at 

(aq, Vo)t - (Vu, • v, Qhu)T + {^b ■ n, u)ot = v). 

Adding sriQhQjV) to both sides of the above equation and then using the definition 
of Qh, we have 

ST(Q?iq,v) + (a(Qoq),vo)r - (Vu, • w,Qhu)T + (v& • n,u)dT = + sriQhfl,^)- 


Now subtracting ( |3.5[ ) from the above equation gives ( |5.7| ). This completes the proof 
of the lemma. □ 

Let T G 77i be an element with e as an edge/face. For any function ip G H^{T), the 
following trace inequality has been derived for general polyhedral partitions satisfying 
the shape regularity assumptions A1 - A4 (see m for details): 

(5.8) \\p\\l < {hf}\\p\\l, + hTW^pWr) ■ 


Theorem 5.3. Under the assumptions of Theorem \ 5.1\ we have the following 
error estimate: 


(5.9) 


( 5] hrWQbU - A.lUr ) < (lklU +2 + llq|U+i). 

XTeTfe / 


Proof By letting v = {0,(5ft,n}, we have vq = 0 and {Vw^r,eh)T = {5h,eh)dT- I* 
follows from (5.7) that 

■ n)gT 

= -Srieh, v) + (Vu, • V, €h)T + sriQhH, v) + £t(m, v) 

= /^^^(eo ' ^h)Q'r T (^/l: ^h)QX 

-hriQotl ■ n - Qb{(l ■ n), Sh)gT “ QhU - u)qt- 
Thus, from the usual Cauchy-Schwarz inequality, we obtain 
, f . < hrlleo ■ n - ef n||gT + ||e/i||aT 

+ ^TllQoq • n - Qbiq ■ n)||aT + \\QhU - u\\ot. 


By first applying the trace inequality (5.8) to the last three terms on the right-hand 
side of ( |5.10[ ), and then summing over all the element T G 7)i we obtain 

^TWShWlr^h^ ^T||eo-n-eb-n|||.r+ XI + ^tIIV e/,||r) 

TeTh TeTh TeTh 

IQoq ~ qllr + ^Tll'^(f 3 oq ~ q)llT) 


(5.11) 


Ten 

+ 'y ' ~ ^IIt T ^tII V (fihU — '»)||r) 


Ten 

< hmebf - 


\£hf + h^W^hWih + 


2 

fc+1 


lfc-l-2 


)■ 


Combining the above inequality with the error estimates 
This completes the proof of the theorem. □ 


and (5.4) yields (5.9). 
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6. Numerical Experiments. In this section, we present some numerical results 
for the hybridized WG-MFEM scheme (3.5 l-( [3?7j) based on the lowest order element 
(i.e., k = 0) for the second order elliptic problem (1.11-( 1.2). Recall that for the lowest 
order hybridized WG-MFEM, the corresponding finite element spaces are given by 

Mh= (g) Vo(T) X Wi(T), 

Ten 

Ah = ■ Mie G Poie),e G £h}- 

For any given v = {vo,v;,}, the discrete weak divergence • v S Pi{T) is 
computed locally on each element T as follows: Find • v G Pi(T) such that 

(Vu, • v,t)t = -(vo, Vt)t -f (Vb ■ n,T)gj,, Vr S Pi(T). 

Let (q/i, Uh, ^h) G A4h X Ah be the hybridized WG-MFEM approximate solution 
arising from (3.5 )-( [3?7| ) and (q, u) be the exact solution of (1.1|-(1.2), respectively. In 
our numerical experiments, we compare the numerical solutions with the projection 
of the exact solution in various norms for their difference: 

~ {^0; ~ Q/iQ Q/i; ^h ~ ^h: ~ Qb'^ ^h- 

In particular, the following norms are used to measure the scale of the error: 


TL^-norm: \\eh\\i,h = ( l|Ve;j|| 
XTen 


^T + h-^ 


E 


L^-norm: 


l|e.|||= E 


\eo\‘^dx -I- hr 


|eo • n — ef, • npds 


' Ten' 


IdT 


L^-norm: \\6h\\ = f E ^^tUh - Qbu\\lT^gfg 


L^-norm: 


Fh = 


Ten 

(s 

^ Ten • 


Ichfdx 


6.1. Example 1. In this test case, the domain is the unit square 12 = (0,1) x 
(0,1) and the coefficient is a = The exact solution is given by u = 

sin(7rx) sin(7rj/). 

Table 6.1 

Example 1. Convergence rate on triangular elements. 


h 

ll|e;i||| 

order 


order 


order 

Ik.ll 

order 

1/4 

2.93e-l 

- 

3.10e-2 

- 

1.22 

- 

1.91e-l 

- 

1/8 

1.47e-l 

0.99 

8.17e-3 

1.92 

5.07e-l 

1.27 

4.76e-2 

2.01 

1/16 

7.39e-2 

1.00 

2.07e-3 

1.98 

2.39e-l 

1.08 

1.19e-2 

2.00 

1/32 

3.70e-2 

1.00 

5.20e-4 

1.99 

1.18e-l 

1.02 

2.97e-3 

2.00 

1/64 

1.85e-2 

1.00 

1.30e-4 

2.00 

5.87e-2 

1.01 

7.42e-4 

2.00 

1/128 

9.25e-3 

1.00 

3.25e-5 

2.00 

2.93e-2 

1.00 

1.86e-4 

2.00 
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This numerical experiment was performed on uniform triangular partitions for 
the domain. The triangular partitions are constructed as follows: 1) first uniformly 
partition the domain into nxn sub-rectangles; 2) then divide each rectangular element 
by the diagonal line with a negative slope. The mesh size is denoted by /i = 1 /n. Table 
|6.1| shows the convergence rate for the hybridized WG-MFEM solutions measured in 
different norms. The results show that the hybridized WG-MFEM approximates are 
convergent with rate 0{h) in and 0{h?) in norms. 

Recall that the primal variable u is approximated by the Lagrange multiplier \h 
as a piecewise constant function on the wired-basket £h- The error function J/i, which 
is the difference of Xu and the projection of the exact solution on each edge e G £h, 
is shown to be convergent at the rate of 0{h^). Thus, the Lagrange multiplier is a 
superconvergent approximation of the exact solution at the mid-point of each edge 
e G £h- This result is similar to the hybridized mixed finite element method. 


6.2. Example 2. In our second test, the domain is again the unit square. 
But the finite element partitions are given by the uniform rectangular meshes. The 
coefficient a is given by a = 1. The data is chosen so that the exact solution is 


u = sin(7ra:) cos(7r?/). The numerical results are presented in Table 6.2 


Table 6.2 

Example 2. Convergence rate on rectangular elements. 


h 

ll|e/i||| 

order 


order 

hh\\i,h 

order 


order 

1/4 

9.19e-l 

- 

1.86e-2 

- 

2.08 

- 

8.70e-l 

- 

1/8 

4.84e-l 

0.93 

4.57e-3 

2.03 

2.09 

- 

2.85e-l 

1.61 

1/16 

2.45e-l 

0.98 

1.14e-3 

2.01 

1.33 

0.65 

7.90e-2 

1.85 

1/32 

1.23e-l 

1.00 

2.84e-4 

2.00 

7.29e-l 

0.87 

2.06e-2 

1.94 

1/64 

6.15e-2 

1.00 

7.10e-5 

2.00 

3.78e-l 

0.95 

5.25e-3 

1.97 

1/128 

3.08e-2 

1.00 

1.77e-5 

2.00 

1.92e-l 

0.98 

1.31e-3 

2.00 


Table |6.2| shows the optimal rate of convergence for the numerical solution in 
and norms. Once again, we see a superconvergence for the solution on the 
wired-basket. 


6.3. Example 3. The test problem here is the same as Example 2, but the 
finite element partitions consist of quadrilaterals constructed as follows. Starting 
with a coarse quadrilateral mesh shown as in Figure [6T] (Left), we successively refine 
each quadrilateral by connecting its barycenter with the middle points of its edges, 
shown as in Figure 6.1 (Right). The numerical results are presented in Table 6.3 All 


the numerical results are in consistency with the theory developed in this paper. 
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Fig. 6.1. Mesh level 1 (Left) and mesh level 2 (Right) for example 3. 


Table 6.3 

Example 3. Convergence rate on quadrilateral elements. 


h 

ll|e/i||| 

order 

II4II 

order 

hh\\l,h 

order 

Ik.ll 

order 

2.86e-l 

8.89e-l 

- 

3.10e-2 

- 

2.03 

- 

8.08e-l 

- 

1.43e-l 

4.74e-l 

0.91 

6.30e-3 

2.30 

1.95 

0.56 

2.71e-l 

1.58 

7.16e-2 

2.42e-l 

0.97 

1.44e-3 

2.13 

1.27 

0.62 

7.66e-2 

1.82 

3.58e-2 

1.21e-l 

1.00 

3.51e-4 

2.04 

7.08e-l 

0.84 

2.03e-2 

1.92 

1.79e-2 

6.03e-2 

1.01 

8.71e-5 

2.01 

3.71e-l 

0.93 

5.25e-3 

1.95 

8.95e-3 

2.95e-2 

1.03 

2.17e-5 

2.00 

1.90e-l 

0.97 

1.36e-3 

1.95 
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